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ABSTRACT 


\l 

To  understand  the  generation  and  the  nonlinear  evolution  of  the  electrostatic 

waves  observed  during  the  AMPTE  (Active  Magnetosphere  Particle  Tracer  Explorers) 

-  ■ - - — — ^ 

solar  wind  releases  e  detaitei^  investigation  is  conducted.  Previous  linear  studies 
ijave^ndicated^tfcat^two  distinct  sets  of  instabilities  may  be  responsible  for  the  ^ 
generations  thoco  woveg.  One  set  consists  of  ion-acoustic  type  instabilities  which 
are  insensitive  to  the  presence  of  a  background  magnetic  field,  while  the  other  group 
corresponds  to  the  modified  two  stream  instabilities  and  requires  the  solar  wind 
flow  to  be  across  the  ambient  magnetic  field,  in-ord^r  to  establish  which  set  of 
instabilitiesAar^more  viably  for  the  generation  of  the  observed  electrostatic  waves 
^a  detailed  linear  Vlasov  theory  has  been  conducted  by  numerically  solving  the  full 
electromagnetic  dispersion  relation.  In  addition  both  the  plasma  wave  as  well  as  the 
magnetic  field  measurements  by  the  IRM  (Ion  Release  Module)  spacecraft  were  used 
to  correlate  the  frequency  and  the  power  of  the  observed  waves  with  the  magnitude 
and  the  direction  of  the  solar  wind  magnetic  field.  The  results  of  these  analyses 
indicate  that  the  ion-acoustic  type  instabilities  have  growth  rates  that  are  an  order 
of  magnitude  or  more  larger  than  those  of  the  modified  two  stream  instabilities. 
Similarly  the  observations  show  no  correlation  between  the  magnitude  and  the 
direction  of  the  ambient  magnetic  field  on  the  one  hand,  and  the  wave  frequency  and 
amplitude  on  the  other  hand.  thus,  indicating  that  the  ion-acoustic  type  instabilities 
are  the  likely  generation  mechanism.  In  order  to  investigate  the  nonlinear  evolution 

y 

of  these  instabilities,  and  discern  their  role  in  the  coupling  of  the  released  ions^to^ 
the  solar  wind,  one-dimensional  electrostatic  simulations  have  been  performed.  The 
results  show  that  both  the  solar  wind  protons  and  the  released  ions  may  be  heated 
and  accelerated  in  the  directions  oblique  to  the  solar  wind  flow  velocity. 


1.  INTRODUCTION 


One  of  the  interesting  and  fundamental  problems  in  the  field  of  space  plasma 
physics  is  the  collisionless  interaction  between  a  flowing  plasma  and  newly  created 
ions  due  to  ionization  of  neutral  gases.  This  type  of  a  problem  can  be  present  under 
a  variety  of  physical  settings  such  as  during  the  interaction  between  the  solar  wind 
and  comets  or  interstellar  helium  gas.  Similarly,  there  have  been  a  number  of  active 
experiments  in  space  in  which  understanding  of  the  interaction  between  the  back¬ 
ground  plasma  and  the  newly  ionized  gases  has  been  a  primary  concern.  Among  these 
are  the  experiments  conducted  during  the  AMPTE  (Active  Magnetospheric  Particle 
Tracer  Explorers)  mission  in  which  lithium  or  barium  gasses  were  released  in  both 
the  solar  wind  and  the  magnetosphere.  One  of  the  main  objectives  of  these  experi¬ 
ments  was  to  trace  the  motion  of  the  released  ions  and  thereby  learn  more  about  the 
transport  properties  of  the  solar  wind  and  the  magnetosphere.  Clearly,  since  both  the 
released  ions  and  the  background  plasma  were  collisionless.  any  interaction  between 
the  two  must  have  taken  place  through  the  electromagnetic  forces  which  were  pre¬ 
sent  in  the  medium  prior  to  the  release  as  well  as  those  that  were  generated  after 
the  release  due  to  collective  plasma  processes.  Thus,  in  order  to  gain  an  understan¬ 
ding  of  the  particle  motion,  it  is  necessary  to  investigate  the  role  of  collective 
plasma  processes.  In  this  study,  the  role  of  wave  particle  interactions  during  the 
AMPTE  solar  wind  releases  is  investigated. 

During  the  AMPTE  mission,  four  releases  were  made  in  the  solar  wind  .  Of 
these  four,  two  were  lithium  releases  which  took  place  on  September  1 1  and  20. 
1984;  and  two  were  barium  releases  which  took  place  on  December  27.  1984,  and 
July  18.  1985.  During  all  of  these  releases,  intense  bursts  of  electrostatic  noise 
were  observed  on  the  upstream  side  of  the  ion  cloud  outside  of  the  diamagnetic 
cavity  [Gurnett  et  al..  1985;  Hausler  et  al..  1986;  Gurnett  et  al.,  1986a.b].  These 
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waves  ranged  in  frequency  between  tens  of  Hz  to  several  KHz  and  their  peak  inte¬ 
grated  electric  field  strength  was  as  high  as  40  mV/m.  An  example  of  these  waves 
is  given  in  the  bottom  panel  of  Figure  1.  In  this  panel  the  electric  field  intensities 
measured  by  the  plasma  wave  instrument  on  board  the  AMPTE/IRM  (Ion  Release  Module) 
spacecraft  during  the  September  1 1  lithium  release  are  shown.  In  this  spectrogram, 
the  axis  of  ordinates  corresponds  to  16  frequency  channels  ranging  from  31.1  Hz  to 
178  KHz.  while  the  axis  of  abscissas  shows  Universal  Time  (UT).  The  top  three 
panels  in  Figure  1  show  the  magnetic  field  measurements  in  geocentric  solar  ecliptic 
(GSE)  coordinates  during  the  same  time  period  [Luhr  et  al.,  1986].  As  can  be  seen  in 
the  third  panel  the  magnetic  field  strength  falls  from  its  solar  wind  value  (~6nT)  to 
nearly  zero  at  07:25:04  UT  indicating  the  entrance  of  the  IRM  spacecraft  into  the 
diamagnetic  cavity  formed  by  the  lithium  cloud.  As  this  cloud  passes  the  spacecraft, 
the  magnetic  field  begins  to  rise  at  about  07:25:15  UT  reaching  values  as  high  as 
23nT.  This  region  of  enhanced  magnetic  field  corresponds  to  a  pile  up  region 
upstream  of  the  lithium  cloud  (for  a  more  detailed  description  of  the  data  see 
Gurnett  et  al..  1986a  and  Luhr  et  al.,  1986).  The  onset  of  intense  electrostatic 
waves  can  be  seen  in  the  bottom  panel  at  about  07:25:30  UT  reaching  peak  power  at 
07:25:35  UT.  As  has  been  shown  by  Gurnett  et  al.  (1985,  1986a)  and  illustrated  in 
more  detail  later  in  this  paper  the  spectrum  of  these  waves  peaks  at  frequencies 
between  a  hundred  and  a  few  hundred  Hertz. 

In  order  to  understand  the  generation  mechanism  of  the  observed  electrostatic 
noise  a  number  of  linear  studies  have  been  performed  [Gurnett  et  al..  1986a.b:  Brinca 
et  al..  1986;  Papadopoulos  et  al..  1987;  Ma  et  al..  1987].  In  Gurnett  et  al.  (1986a.b) 
it  was  suggested  that  the  wave  amplification  takes  place  via  an  ion-acoustic  type 
instability  due  to  the  relative  drift  between  the  solar  wind  and  the  newly  created 
ions.  In  these  studies  growth  rates  of  the  waves  with  phase  velocity  along  the  solar 
wind  flow  were  calculated  and  shown  to  be  significant.  These  studies  were  later 
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extended  to  obliquely  propagating  waves  showing  that  the  growth  rates  of  these  waves 
are  comparable  or  larger  than  the  parallel  propagating  waves  [Brinca  et  al..  1986;  Ma 
et  ah.  1987].  Specifically,  it  was  shown  by  Ha  et  ah  (1987)  that  two  types  of 
instabilities  are  possible.  One  is  an  instability  due  to  the  relative  drift  between 
the  released  ions  and  all  the  electrons,  while  the  other  is  due  to  the  streaming 
between  the  released  ions  and  the  solar  wind  protons.  Both  of  these  instabilities, 
however,  are  ion-acoustic  type  in  that  the  wave  dispersion  properties  are  similar  to 
that  of  ion-acoustic  waves.  In  fact  using  the  nomenclature  used  by  Gary  and  Omidi 
(1987),  the  first  type  of  instability  corresponds  to  the  electron/ion-acoustic,  while 
the  second  one  corresponds  to  the  ion/ion- acoustic  instabilities.  Of  primary 
importance  is  the  fact  that  these  two  instabilities  are  insensitive  to  the  presence 
of  a  background  magnetic  field  [Akimoto  and  Omidi,  1986;  Gary  and  Omidi,  1987]. 

In  addition  to  the  above  set  of  instabilities,  another  set  has  also  been  suggested 
by  Papadopoulos  et  al.  (1987)  to  be  responsible  for  the  generation  of  the  electro¬ 
static  waves.  These  types  of  instabilities  arise  when  a  relative  drift  between  the 
electrons  and  ions  or  two  ion  species  exists  across  a  magnetic  field.  In  these  insta¬ 
bilities.  it  is  essential  that  the  electrons  be  treated  as  magnetized  while  the  ions 
may  be  considered  as  unmagnetized.  Because  the  magnetic  field  in  a  plasma  plays  an 
important  role  for  these  modes,  the  parameter  0  (the  ratio  of  kinetic  to  magnetic 
pressure)  can  greatly  affect  the  growth  rates  of  these  waves.  In  fact  while  these 
waves  are  adequately  described  by  the  electrostatic  dispersion  relation  at  low  values 
of  0,  electromagnetic  effects  become  important  as  0  increases,  until  eventually  an 
electromagnetic  dispersion  relation  must  be  used  to  discern  their  properties  such  as 
growth  rates. 

Although  in  general  it  is  possible  for  all  of  the  above  mentioned  instabilities  to 
be  operative  simultaneously,  it  is  much  more  likely  that  one  of  them  becomes  domi¬ 
nant  and  thereby  governs  the  nature  of  the  wave-particle  interactions  during  the  non- 


tinear  regime.  Because  our  main  objective  is  to  discern  the  role  of  the  observed 
waves  on  particle  dynamics,  it  is  essential  that  the  likely  instability  responsible  for 
the  generation  of  these  waves  be  identified.  This  is  the  aim  of  the  next  section  in 
which  the  growth  rate  calculations  of  Papadopoulos  et  al.  (1987)  will  be  extended  by 
using  the  full  electromagnetic  dispersion  relation.  In  addition,  the  plasma  wave  and 
the  magnetic  field  measurements  will  be  utilized  to  discern  the  role  of  the  magnetic 
field  magnitude  and  direction  on  the  wave  properties,  and  thereby  deduce  the  nature 
of  the  instability.  As  will  be  shown  in  the  next  section  the  electron/ion-acoustic 
and  the  ion/ion-acoustic  instabilities  were  more  likely  responsible  for  the  generation 
of  the  electrostatic  waves. 

Having  determined  the  nature  of  the  instability  in  Section  2.  a  detailed  simula¬ 
tion  study  will  be  conducted  in  Section  3  to  discern  the  nonlinear  evolution  of  the 
waves  and  thereby  gain  a  better  understanding  of  the  wave-particle  interactions 
during  the  AMPTE  solar  wind  releases.  The  simulations  consist  of  two  parts:  one  in 
which  the  electrons  and  the  ions  are  treated  as  particles  with  unrealistic  mass 
ratios.  The  results  of  such  a  simulation  study  were  briefly  described  by  Omidi  et 
al.  (1987).  The  other  type  of  simulation  consists  of  fluid  electrons  and  particle 
ions.  The  advantage  of  this  kind  of  simulation  method  is  that  realistic  mass  ratios 
are  used  and  thereby  the  observed  plasma  parameters  can  be  utilized.  A  summary  and 
discussion  of  the  results  will  be  presented  in  Section  4. 


2.  LINEAR  ANALYSIS 


2.1  Theoru 

As  mentioned  earlier  one  way  in  which  the  viability  of  the  various  proposed 
instabilities  in  generating  the  observed  electrostatic  noise  can  be  tested  is  by  com¬ 
paring  their  linear  growth  rates.  To  do  so  we  use  a  plasma  model  similar  to  that 
given  by  Gurnett  et  al.  (1986a)  in  which  the  released  ions  (in  this  case  lithium) 
form  a  cold  non-drifting  Maxwellian  with  a  temperature  TLi  =  2  *  103°K.  The  solar 
wind  protons  form  a  Maxwellian  with  temperature  Tp  =  1 0S°K  and  a  drift  speed 
Vsw  =  450  km/s.  As  for  the  electrons  two  possibilities  were  allowed;  in  one  case 

the  electrons  consisted  of  two  separate  Maxwellian  populations  with  temperatures 
h  c 

Te  =  5*1 0S°K  and  Te  =  5*104°K  corresponding  to  solar  wind  electrons  and  those  asso¬ 
ciated  with  the  photoionization  of  the  released  gas  respectively.  In  the  other  case 
only  the  solar  wind  electrons  were  present  allowing  for  the  fact  that  the  colder 

photoelectrons  get  replaced  by  the  hotter  solar  wind  electrons.  In  both  cases  the 
electrons  were  assumed  to  have  a  drift  speed  such  that  the  net  current  was  zero.  It 
should,  however,  be  mentioned  that  growth  rate  calculations  with  finite  current  in 
the  plasma  were  also  performed  showing  no  significant  difference  in  the  results  (see 
Gurnett  et  al..  1986a).  Finally  in  regards  to  the  density  ratio  of  the  released  ions 
to  the  solar  wind  protons  a  wide  range  of  values  were  considered.  Observations  show 
that  the  density  ratio  changes  considerably  with  distance  from  the  outer  edge  of  the 
diamagnetic  cavity.  The  results  of  linear  theory  have  shown  that  the  growth  rates 
are  maximum  when  the  proton  and  the  lithium  (or  barium)  densities  are  equal. 
Throughout  this  paper,  we  take  the  density  ratio  to  be  one,  because  the  role  of  this 
parameter  is  the  same  in  both  sets  of  instabilities. 

Because  the  growth  rate  calculations  of  the  ion-acoustic  type  instabilities  have 
been  performed  in  great  detail  in  previous  studies  [Gurnett  et  al.,  1986a,b;  Brinca  et 


al.,1986;  Ma  et  al.,  1987],  here  we  will  only  summarize  the  pertinent  results. 

Figure  2  (from  Ma  et  al..  1987)  shows  the  growth  rates  (2T max)  of  electrostatic 
waves  maximized  over  the  wave  number  «)  as  a  function  of  9  (the  angle  between  the 
K  vector  and  the  solar  wind  velocity)  for  the  parameters  shown  in  the  figure.  The 
solid  line  represents  the  case  where  no  cold  photoelectrons  are  present,  while  the 
dashed  line  corresponds  to  the  case  where  photoelectrons  are  present.  As  can 
be  seen,  the  growth  rates  decrease  monotonically  with  increasing  e  for  e  <  70°. 
However,  beyond  this  point  and  in  the  range  70°  <  9  <  82°  the  growth  rates  rise  with 
9.  indicating  a  change  in  the  nature  of  the  instability.  As  was  pointed  out  by  Ma  et 
al.  (1987).  the  growth  rates  at  9  <  70°  correspond  to  a  kinetic  instability  due  to  the 
relative  drift  between  the  electrons  and  the  lithium  ions.  This  instability  is  simi¬ 
lar  to  a  current  driven  ion-acoustic  instability  in  that  the  wave  growth  is  due  to 
Landau  resonance  with  the  electrons.  Using  the  terminology  introduced  by  Gary  and 
Omidi  (1987).  we  refer  to  this  instability  as  the  electron/ion-acoustic  instability. 
The  growth  rates  at  9>  70°.  on  the  other  hand  correspond  to  a  fluid  type  instability 

(for  the  temperatures  considered  here)  due  to  the  relative  drift  between  the  two  ion 
species.  This  instability  was  referred  to  as  the  ion/ ion- acoustic  instability  by  Gary 

and  Omidi  (1987). 

It  is  evident  from  Figure  2  that  in  the  absence  of  photoelectrons,  the  ion/ion- 
acoustic  instability  has  a  maximum  growth  rate  which  is  larger  than  that  of  the 
electron/ion-acoustic  by  an  order  of  magnitude.  Thus,  in  the  absence  of  cold  elec¬ 
trons  the  ion/ion-acoustic  instability  is  probably  the  dominant  one.  On  the  other 
hand,  the  growth  rates  of  the  two  instabilities  are  comparable  when  both  electron 
species  are  present.  This  effect  comes  about  by  virtue  of  the  fact  that  the  electron/ 
ion-acoustic  waves  grow  due  to  Landau  resonance  with  the  electrons  whose  distri¬ 
bution  function  attains  a  larger  slope  in  the  presence  of  the  cold  electrons. 

Similarly  the  ion/ion-acoustic  waves  suffer  a  greater  Landau  damping  due  to  the 


increase  in  the  slope  of  the  electron  distribution  function.  As  was  mentioned  earlier 
these  two  instabilities  are  insensitive  to  the  presence  of  a  background  magnetic 
field,  and  thus  one  would  not  expect  to  see  any  correlation  between  the  electrostatic 
waves  observed  during  the  AMPTE  release  and  the  magnetic  fluctuations  in  the  solar 
wind. 

The  other  set  of  instabilities  which  have  been  proposed  to  explain  the  electro¬ 
static  noise  are  the  modified  two  stream  instabilities  [Papadopoulos  et  al.,  1987]. 
These  instabilities  also  arise  due  to  either  a  drift  between  the  electrons  and  the 
ions  or  two  ion  species.  However,  in  either  case  the  drift  must  be  oblique  to  the 
direction  of  the  magnetic  field  with  the  growth  rates  becoming  maximum  when  the 
drift  is  perpendicular.  To  help  better  visualize  the  geometry  of  these  instabilities 
we  show  in  Figure  3  the  coordinate  system  (p,  t\,  C  )  in  which  the  solar  wind  flow 
velocity  (Vsw)  is  along  the  p-axis  and  the  ambient  magnetic  field  (B0)  is  in  the  C 
direction.  Also  shown  in  Figure  3  is  a  typical  wave  vector  K  which  makes  an  angle  <t> 
with  the  C-axis  and  its  projection  on  the  (p.T\)  plane  makes  an  angle  e  with  the 
p-axis.  An  important  point  to  note  is  that  unlike  the  ion-acoustic  instabilities,  the 
modified  two  stream  instabilities  are  not  axially  symmetric  with  respect  to  the 
p-axis  and  thus  both  angles  9  and  $  play  a  role  in  the  growth  rates  of  the  waves. 
Another  important  property  of  these  instabilities  is  the  fact  that  while  in  the  low 
0  regime  they  are  well  described  by  the  electrostatic  approximation,  at  larger  values 
of  0  the  electromagnetic  effects  become  important  and  can  lead  to  their  stabiliza¬ 
tion  (for  more  details  on  the  modified  two  stream  instabilities  see,  e.g.  Gary  et  al., 
1987).  Since  in  the  analysis  of  Papadopoulos  et  al.  (1987)  an  approximate  dispersion 
relation  was  utilized  and  also  the  value  of  0  was  taken  to  be  zero,  it  is  essential 
that  their  study  be  extended  by  using  a  full  electromagnetic  dispersion  relation  and 

allow  for  finite  values  of  0. 

The  general  electromagnetic  form  of  the  linear  Vlasov  dispersion  equation  for  a 
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iwo-dimensional  wave  vectorOC)  is  given  in  Tokar  and  Gary  (1985).  Using  this  dis¬ 
persion  equation  and  the  plasma  parameters  described  earlier,  growth  rates  have 
been  calculated  for  a  wide  range  of  electron  0  (0e)  and  angles  e  and  <t>.  The  results 
of  these  computations  for  the  same  parameters  shown  in  Figure  2  but  without  the 

cold  electrons  are  shown  in  Figures  4-6.  In  Figure  4  the  growth  rates  normalized  to 
the  proton  angular  gyrofrequency  (Qj)  and  maximized  over  K  are  shown  as  a  function 
of  4>  for  a  number  of  £e  values  at  9  =  0°.  As  was  mentioned  earlier,  two  types  of 
instabilities  are  possible  when  the  flow  is  across  the  magnetic  field.  One  is  due  to 
the  relative  motion  between  the  electrons  and  ions  and  is  referred  to  as  the 
electron/ion  cross-field  instability  by  Gary  et  al.  (1987),  while  the  other  is  due  to 
the  relative  drift  between  the  ions  and  is  referred  to  as  the  ion/ion  cross-field 

instability.  The  growth  rates  shown  at  $  <  88°  are  due  to  the  electron/ion  cross¬ 
field  instability  while  those  at  <!>  >  88°  are  due  to  the  ion/ion  cross-field  instabi¬ 
lity.  As  can  be  seen  in  Figure  4.  the  latter  mode  becomes  stable  as  £e  is  increased 
to  about  0e  -0.1.  It  is  also  evident  that  as  fie  increases  the  maximum  growth  rates 
decrease  indicating  that  the  electromagnetic  effects  become  important  and  tend  to 
stabilize  both  waves.  It  should,  however,  be  mentioned  that  while  the  ion/ion  mode 
becomes  stable  at  Q  =  0  when  0e  becomes  large,  it  still  remains  unstable  at  other 
values  of  0.  To  demonstrate  this  point,  we  show  in  Figure  (5a)  the  maximum  growth 
rates  as  a  function  of  0  at  $  s  88.5°  for  four  different  values  of  0e.  The  growth 
rates  shown  in  this  panel  correspond  to  the  ion/ion  cross-field  instability.  As  can 
be  seen,  when  £e  is  small  (0.01)  maximum  growth  rates  peak  at  0  =  0°.  However,  as 
0e  increases  the  peak  in  2fmax  shifts  to  higher  values  of  0.  Note  that  this  effect  is 
in  contrast  to  that  of  the  electron/ion  mode  whose  growth  rates  are  shown  in  panel 
(b)  of  Figure  5.  As  can  be  seen  in  this  panel  the  growth  rates  drop  monotonically 
with  0  regardless  of  the  value  of  0e.  Again  it  is  evident  in  panel  (a)  of  Figure  5 
that  the  maximum  growth  rate  of  the  ion/ion  mode  decreases  with  increasing  £e  due 
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to  the  electromagnetic  effects. 

A  point  that  should  be  mentioned  in  connection  with  Figure  (5b)  is  that  the 
growth  rates  at  0e  =  0.1  are  larger  than  those  at  0e  =  0.5  and  in  turn  those  growth 
rates  are  larger  than  the  ones  at  0e  =  0-01-  This  effect  is  due  to  the  fact  that  in 
Figure  (5b)  $  is  kept  fixed  while  3e  has  changed.  However,  as  can  be  seen  in  Figure 
4,  the  value  of  <t>  at  which  2fmax  peaks  varies  with  3e.  In  short  contrary  to  what 
might  appear  in  panel  (b)  of  Figure  5.  the  maximum  growth  rates  of  both  of  the 
instabilities  decrease  with  increasing  0e. 

Thus  far  the  general  behavior  of  the  two  cross-field  instabilities  have  been 
investigated  for  various  values  of  3e.  ln  order  to  determine  the  viability  of  these 
instabilities  in  generating  the  observed  electrostatic  noise,  it  is  best  to  look  at  the 
real  and  the  imaginary  parts  of  the  frequency  as  a  function  of  K  for  a  representative 
value  of  0e.  Using  the  plasma  parameters  given  earlier  (ne=  16cm-i5.  Te=  5  x  10s*  K) 
and  a  magnetic  field  strength  of  23  (nT)  (see  Figure  1 )  we  obtain  the  value  of  0e  -  0.5. 
Note  that  this  value  is  probably  a  minimum  in  that  the  magnetic  field  strength  used 
in  its  computation  is  the  maximum  observed  value.  Using  this  value  of  0e.  the  real 
(<or)  and  the  imaginary  (2T)  parts  of  the  frequency  have  been  calculated  as  a  function 
of  K  for  values  of  e  and  $  at  which  maximum  growth  occurs.  Note  that  the  compu¬ 
tations  were  carried  out  in  a  frame  where  the  lithium  ions  are  at  rest  which  should 
correspond  to  that  of  the  satellite.  Panel  (a)  of  Figure  6  shows  o>r  (solid)  and  % 
(dashed)  for  the  ion/ion  cross-field  instability  while  panel  (b)  shows  those  of  the 
electron/ion  cross-field  instability  as  a  function  of  the  wave  number  (K)  times  the 
proton  gyroradius  (p^).  Using  the  maximum  magnetic  field  strength  observed  (23  nT) 
during  the  September  1 1  release,  it  is  straightforward  to  compute  the  maximum 
growth  rates  of  the  two  instabilities  and  their  corresponding  frequency  and  wave¬ 
length.  As  can  be  seen  in  Figure  6  the  maximum  growth  rate  of  the  ion/ion  mode  is 
V  ~  4Qj  -  9  Hz,  while  that  of  the  electron/ion  mode  is  Z  =  5Qi  =  1 1 .3  Hz.  A  compar- 
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ison  between  these  growth  rates  and  those  shown  in  Figure  2  clearly  demonstrate 
that  the  growth  rates  of  the  ion-acoustic  instabilities  are  more  than  an  order  of 
magnitude  larger  than  those  of  the  cross-field  modes.  It  is  also  evident  from  Figure 
6  that  the  maximum  growth  of  the  ion/ion  cross-field  instability  occurs  at  a  fre¬ 
quency  f  =  1 7.5  Hz  while  that  of  the  electron/  ion  occurs  at  f  =  52  Hz.  Note,  however, 
that  these  frequencies  are  too  low  to  explain  the  observed  spectrum.  Finally  both  of 
the  instabilities  have  maximum  growth  at  a  wave  length  of  about  X  =  100  km.  These 
wavelengths  also  seem  to  be  too  large  in  view  of  the  fact  that  they  are  comparable 
or  larger  than  the  entire  size  of  the  lithium  cloud.  Thus  the  results  of  linear  theory 
indicate  that  the  ion-acoustic  instabilities  are  much  more  likely  to  have  been 
responsible  for  the  generation  of  the  observed  electrostatic  noise. 

2.2  Observations 

In  addition  to  linear  theory,  the  measurements  by  the  plasma  wave  instrument 
and  the  magnetometer  can  also  be  utilized  to  assess  the  plausibility  of  each  set  of 
instabilities  having  resulted  in  the  generation  of  the  electrostatic  waves.  This  can 
be  done  by  noting  that  while  the  ion-acoustic  type  instabilities  are  insensitive  to  the 
presence  of  an  ambient  magnetic  field,  the  cross-field  instabilities  result  in  wave 
frequencies  and  growth  rates  which  are  sensitive  to  both  the  magnitude  and  the  direc¬ 
tion  of  the  magnetic  field.  Specifically  one  would  expect  the  wave  frequencies  of 
the  cross-field  instabilities  to  increase  with  the  ambient  magnetic  field  strength, 
and  also  the  wave  power  to  rise  (fall)  as  the  magnetic  field  direction  becomes  more 
(less)  oblique  to  the  solar  wind  flow  direction.  Note  also  that  one  would  not  expect 
the  cross-field  instabilities  to  be  operative  when  the  solar  wind  flow  is  parallel  or 
quasi-parallel  to  the  magnetic  field.  The  September  1 1th  release  is  particularly 
suitable  for  conducting  this  test  in  that  both  the  magnitude  and  the  direction  of  the 
magnetic  field  changed  considerably  during  the  period  when  the  electrostatic  waves 


were  observed  (Lii'hr  et  al.,  1986).  For  example,  it  is  evident  in  the  third  panel  of 
Figure  1  that  the  magnetic  field  varies  from  B  =  8nT  at  07:25:30  UT  to  B  =  23nT  at 
07:25:35  UT  and  then  drops  back  to  B-IInT  at  07:25:38  UT.  Similarly,  assuming  the 
solar  wind  flow  to  be  along  the  Sun-Earth  line  in  the  GSE  coordinate  system  (i.e.  4>B 
=  360°  and  e0  =  0a),  it  is  clear  from  the  top  two  panels  in  Figure  1  that  the  mag¬ 
netic  field  direction  changes  considerably  with  respect  to  the  flow  velocity. 

The  eight  panels  in  Figure  7  show  the  electric  field  spectrum  densities  during 
the  period  07:25:30.2  -  07:25:38  UT.  The  begining  of  this  period  corresponds  roughly 
to  the  onset  of  the  electrostatic  waves  with  the  power  reaching  its  maximum  at 
about  07:25:35.  Beyond  this  time,  the  wave  amplitudes  decrease  approaching  the 
background  level  at  around  07:25:40  UT.  One  striking  feature  of  the  power  spectrums 
shown  in  Figure  7  is  the  fact  that  while  the  shape  of  the  spectrums  change  with  time, 
the  maximum  power  remains  at  a  relatively  fixed  frequency  of  f-150  Hz.  This 
feature  argues  against  the  cross-field  instabilities  being  responsible  for  the  ampli¬ 
fication  of  the  waves  in  two  ways.  First,  the  frequencies  at  peak  power  are  much 
larger  than  the  frequencies  (f  <50  Hz  assuming  B  =  23nT)  at  which  cross-field 
instabilities  have  maximum  growth  rates.  Second,  since  the  magnetic  field  strength 
varies  considerably  during  the  time  interval  considered  in  Figure  7.  one  would  expect 
the  frequencies  at  peak  power  to  vary  accordingly.  This  effect,  however,  is  not 
reflected  in  the  power  spectrums  shown  in  Figure  7.  In  regards  to  the  direction  of 
the  magnetic  field  it  is  evident  from  the  top  two  panels  of  Figure  1  that  both  $B 
(azimuthal  angle)  and  eB  (latitudinal  angle)  fluctuate  by  as  much  as  30°  on  relatively 
short  time  scales.  These  fluctuations  however,  do  not  seem  to  affect  the  intensity 
of  the  waves.  In  addition,  while  the  magnetic  field  may  at  times  be  oblique  to  the 
flow  direction,  it  is  far  from  being  perpendicular  to  it.  Again  suggesting  that  the 
prevailing  conditions  in  the  plasma  are  not  conducive  to  the  generation  of  waves  via 
the  cross-field  instabilities. 


In  short  both  theoretical  investigations  as  well  as  observational  evidence  indi¬ 
cate  that  the  electrostatic  waves  observed  during  the  AMPTE  solar  wind  releases 
were  generated  by  the  ion-acoustic  type  instabilities. 
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3. 


NONLINEAR  ANALYSIS  ! 


! 


3.1  General  Remarks 

Thus  far  we  have  considered  two  sets  of  instabilities  with  the  main  distinction 
between  the  two  being  that  one  set  requires  the  presence  of  a  background  magnetic 
field  while  the  other  can  exist  in  an  unmagnetized  plasma.  Having  shown  that  the 
ion-acoustic  type  instabilities  are  the  dominant  modes,  the  remainder  of  this  paper 
will  be  devoted  to  the  nonlinear  evolution  of  the  waves.  A  point  to  note  is  that 
while  the  electron/ion-acoustic  and  the  ion/ion-acoustic  instabilities  result  in 
plasma  waves  with  similar  dispersive  properties  (i.e.  acoustic  like),  it  is  nonethe¬ 
less  essential  that  these  two  instabilities  be  distinguished  when  considering  their 
nonlinear  evolution.  This  is  because  the  source  of  free  energy  and  therefore,  the 
saturation  mechanism  of  each  instability  is  different.  Specifically,  the  electron/ 
ion-acoustic  waves  grow  due  to  Landau  resonance  with  the  electrons  and  therefore, 
are  expected  to  saturate  via  quasi-linear  diffusion  of  the  electrons  and  formation  of 
a  plateau  in  their  distribution  function.  Because  the  time  scales  associated  with 
this  process  are  small  as  compared  to  ion  time  scales,  one  would  not  expect  the 
electron/ion-acoustic  instability  to  affect  the  ions  in  an  appreciable  way.  On  the 
other  hand  the  ion/ion-acoustic  waves  grow  due  to  interaction  between  the  two  ion 
species  and  thus  their  saturation  mechanism  must  involve  ions.  This  can  also  be  seen 
by  noting  that  any  heating  of  the  electrons  by  the  waves  will  tend  to  increase  the 
ratio  of  Te/Tp  which  enhances  the  growth  rates  (see  e.g.  Gary  and  Omidi.  1987).  In 
other  words,  electron  heating  does  not  have  a  stabilizing  effect  on  the  ion/ion- 
acoustic  waves. 

As  can  be  seen  in  Figure  2  when  both  the  cold  photoelectrons  as  well  as  the  hot 
solar  wind  electrons  are  present  the  maximum  growth  rates  of  the  electron/ion- 
acoustic  and  the  ion/ion-acoustic  waves  are  comparable.  Thus,  it  is  conceivable  that 
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both  modes  may  grow  with  the  electron/ion  waves  heating  the  cold  electrons  and  the 
ion/ion  waves  interacting  mainly  with  the  two  ion  populations.  On  the  other  hand, 
when  no  cold  electrons  are  present  the  ion/ion  mode  becomes  dominant.  Although  the 
effect  of  waves  on  the  electrons  is  interesting,  our  main  goal  is  to  see  how  the 
newly  born  ions  are  affected  by  the  waves,  that  is  to  say.  how  the  observed  waves 
could  have  resulted  in  coupling  between  the  solar  wind  protons  and  the  lithium  or 
barium  ions.  Thus  we  j.v_  *  only  be  concerned  with  the  ion/ion  mode  and  investigate 
its  nonlinear  evolution.  To  this  end.  one-dimensional  electrostatic  particle  simula¬ 
tions  have  been  performed,  the  results  of  which  are  presented  here.  These  results 
correspond  to  two  classes  of  simulations!  one  in  which  the  electrons  and  the  ions  are 
treated  as  particles  with  the  disadvantage  of  an  unrealistic  mass  ratio  between  the 
electrons  and  the  ions.  The  other  set  of  results  correspond  to  a  simulation  code  in 
which  the  electrons  are  approximated  by  a  Boltzman  fluid.  This  code  has  the  advan¬ 
tage  that  the  electron  to  ion  mass  ratio  is  realistic  and  therefore  the  observed 
plasma  parameters  can  be  used  in  the  simulations. 

?,2  MLEamslfl  Simulations 

The  full  particle  simulation  model  consists  of  a  one-dimensional  periodic  box  in 
the  direction  of  the  K  vector(e.g.  the  x-axis  in  Figure  3  when  $  =  90*).  Note  that 
the  K  vector  can  be  directed  in  any  desired  direction  by  assigning  the  flow  velocity 
of  each  plasma  species  to  be  equal  to  its  original  flow  velocity  projected  along  the 
desired  direction.  Thus  for  example,  in  order  to  simulate  the  waves  propagating  at 
an  angle  e  with  respect  to  the  solar  wind  flow  direction  (Vsw)  it  is  sufficient  to 
assign  the  solar  wind  velocity  to  be  Vswcose.  The  simulation  results  presented  in 
this  paper  all  correspond  to  a  case  where  no  photoelectrons  are  present,  and  thus,  the 
ion/ion  mode  has  growth  rates  that  are  larger  than  those  of  the  electron/ion  mode  by 
an  order  of  magnitude.  This  is  desirable  since  our  primary  interest  is  in  the  ion/ion 


mode.  Also  the  fact  that  the  ion/ion  mode  has  larger  growth  rates  makes  the  use 
of  a  one-dimensional  model  appropriate  in  view  of  the  fact  that  even  in  a  two- 
dimensional  model  the  electron/ion  modes  would  not  have  had  a  chance  to  grow.  This 
point  will  be  substantiated  later  when  it  is  shown  that  the  ions  become  disrupted 
on  time  scales  shorter  than  the  growth  time  of  the  electron/  ion  mode.  The  size  of 
the  simulation  box  is  32  electron  Oebeye  lengths  (Xe).  and  it  consists  of  64  cells. 

The  simulation  frame  is  taken  to  be  the  same  as  that  in  the  linear  theory  and  thus, 
initially  the  lithium  ions  form  a  non-drifting  Maxwellian  while  the  solar  wind 
protons  consist  of  a  Maxwellian  drifting  with  velocity  Vswcose.  Similarly  the 
electrons  have  a  drift  speed  such  that  the  net  current  is  zero. 

As  mentioned  earlier,  one  of  the  difficulties  associated  with  a  full  particle 
code  is  the  fact  that  the  simulation  must  account  for  the  electron  dynamics. 

Therefore,  the  time  steps  in  the  simulations  must  be  small  enough  such  that  the 
electron  dynamics  are  resolved.  On  the  other  hand,  since  the  ion/ion  mode  involves 
the  motion  of  the  ions  which  takes  place  on  a  much  longer  time  scale  it  is  not 
feasible  to  use  a  realistic  mass  ratio  between  the  electrons  and  the  ions.  Thus,  it 
is  necessary  to  adjust  the  ratio  of  the  masses  such  that  both  the  electron  and  the  ion 
dynamics  can  be  simultaneously  resolved.  For  this  reason  we  use  mLj  =  6mp  =  600 
me  where  Li.  p.  and  e  stand  for  lithium,  proton  and  electron  respectively  and  m  is 
the  mass. 

Although  the  use  of  unrealistic  mass  ratios  does  not  usually  change  the  physical 
nature  of  the  problem,  it  does  make  it  impossible  to  use  the  exact  plasma  parameters 
as  those  in  a  realistic  plasma.  The  results  of  the  full  particle  simulation  presented 
here  correspond  to  the  following  plasma  parameters.  The  proton  flow  velocity  Vsw  = 
0.71  ve.  its  thermal  velocity  vp  s  0.01  ve  where  ve  is  the  electron  thermal  speed. 

The  lithium  thermal  velocity  vLi  =  0.26*1 0"2ve,  the  angle  of  propagation  e  =  85*.  and 
the  lithium  density  Ny  =  Np.  In  these  simulations  the  electrons  are  represented  by 
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32.000  particles  while  the  protons  and  the  lithium  ions  are  each  represented  by 
16.000  particles.  The  time  step  At  =  0.2o>pe_1  where  o>pe  is  the  electron  plasma 
frequency.  Finally,  it  should  be  mentioned  that  the  simulations  are  carried  out  as  an 
initial  value  problem  where  the  system  is  allowed  to  evolve  in  time  without  a 
constant  refurbishing  of  new  particles. 

The  four  panels  in  Figure  8  show  the  energy  history  of  the  electrostatic 
waves,  electrons,  protons,  and  the  lithium  ions  from  top  to  bottom  respectively. 

These  energies  are  normalized  to  the  total  initial  energy  in  the  system.  As  can  be 
seen  in  the  top  panel,  initially  the  wave  energy  increases  with  time  till  it  eventu¬ 
ally  saturates  at  about  t~560o>pe-1.  Thus  during  this  period  the  ion/ion  waves  are 
growing  with  time.  Beyond  t  =  SSOOpg"1  the  wave  energy  begins  to  decrease  due  to 
damping  by  wave-particle  interactions.  The  damping  of  the  waves  lasts  until  about  t  - 
700<*)pe_1.  when  the  level  of  fluctuations  becomes  more  steady  with  time.  The 
second  panel  in  Figure  8  shows  the  energy  history  of  the  electrons  which  closely 
mimics  that  of  the  waves.  This  behavior  is  in  contrast  to  the  energy  history  of  the 
protons  and  the  lithium  ions.  As  can  be  seen,  the  protons  which  are  the  source  of 
free  energy  for  wave  amplification,  lose  energy  with  time  until  t  ~  560<i>pe_1. 

Beyond  this  point,  their  energy  remains  relatively  constant  with  time.  The  lithium 
ions  gain  energy  throughout  the  run,  although  the  rate  of  increase  becomes  quite 
small  after  t  =  675o)pe_1.  The  behavior  of  the  protons  and  the  lithium  ions  sub¬ 
stantiate  the  expectation  that  the  ion/ion  mode  can  greatly  affect  the  ions.  One  of 
the  interesting  aspects  of  Figure  8  is  the  fact  that  the  energy  history  of  the  elec¬ 
trons  and  the  waves  are  very  similar.  This  suggests  that  the  observed  changes  in 
the  energy  of  the  electrons  must  be  for  the  most  part  due  to  coherent  motion  in  the 
electric  field  of  the  waves,  rather  than  random  motions  which  cause  an  increase  in 
their  temperature.  This  point  will  be  further  substantiated  later. 

To  better  understand  the  role  of  the  wave-particle  interactions  on  ion  dynamics 
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the  phase  space  density  plots  (Vx  vs  x)  of  the  protons  and  the  lithium  ions  are 
shown  in  Figures  9  and  10  at  three  separate  times  during  the  simulations.  In  these 
figures,  the  velocity  of  each  species  is  normalized  to  its  thermal  speed.  As  can  be 
seen  in  panel  (a)  of  Figure  9,  initially  the  protons  form  a  beam  with  a  drift  speed 


vd  =  0.71cos(85°)  Vq ~  6.2Vp.  At  later  times  when  the  waves  have  grown,  they  begin 
to  interact  with  the  protons  (panel  b)  leading  to  their  eventual  trapping  (panel  c). 
Note  in  panel  (b)  that  six  oscillations  are  present,  which  agrees  with  the  results  of 
linear  theory  indicating  that  modes  5  and  6  are  the  most  unstable  ones.  From  the 
number  of  trapping  vortices  in  panel  (c).  however,  it  is  clear  that  mode  number  5  has 
become  dominant  at  later  times.  Note  also  that  the  center  of  the  trapping  vortices 
lies  around  Vx  -  2vp  which  is  reasonable  in  view  of  the  fact  that  the  phase  velocity 
of  the  unstable  waves  is  V$  =  1.7vp.  The  phase  space  density  plots  of  the  lithium 
ions  presented  in  Figure  10  show  that  initially,  these  ions  have  a  zero  drift  velocity. 
At  later  times  during  the  simulation,  the  lithium  ions  begin  to  interact  with  the 
waves  resulting  in  the  oscillations  shown  in  panel  (b).  Note  again  that  six  oscilla¬ 
tions  can  be  seen  during  this  relatively  early  stage  of  the  simulation.  Panel  (c)  of 
Figure  (10)  shows  the  lithium  ions  becoming  trapped,  indicating  that  the  nonlinear 
saturation  mechanism  of  the  ion/ion  mode  is  in  this  case  trapping  of  both  of  the  ion 
species.  This  behavior  may  at  first  seem  somewhat  surprising  in  that  the  trapping 


frequency  <oT  a  (flit)1  '2  (where  E  is  the  electric  field  amplitude,  q  and  m  are  the 
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charge  and  mass  of  the  particle)  of  each  species  is  proportional  to  the  inverse  square 
root  of  its  mass.  Thus,  ordinarily  one  would  expect  the  protons  to  get  trapped  and 
saturate  the  instability  before  the  lithium  ions  have  had  a  chance  to  get  trapped. 

The  fact  that  in  this  case  both  species  suffered  trapping  can  be  explained  by  noting 
that  the  phase  velocity  of  the  waves  V$  =  1 .7vp,  falls  closer  to  the  lithium  ion  dis¬ 
tribution  as  compared  to  that  of  the  protons  whose  beam  speed  is  Vd  =  6.2vp.  Thus, 
in  order  for  the  protons  to  get  trapped  they  must  first  get  decelerated  by  a  substan- 
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tial  amount,  and  a  large  electric  field  must  be  set  up.  This,  however,  requires  a 
longer  time  which  implies  a  smaller  effective  trapping  frequency  for  the  protons.  In 
other  words  due  to  the  large  flow  velocity  of  the  protons  relative  to  the  phase  speed 
of  the  waves,  their  effective  trapping  frequency  becomes  comparable  to  that  of  the 
lithium  ions  whose  velocity  distribution  function  lies  closer  to  the  phase  velocity  of 
the  waves. 

In  order  to  understand  in  more  detail  the  effect  of  wave-particle  interactions 
on  the  electrons,  their  velocity  distribution  functions  at  t  =  0  and  600u)pe_1  are 
shown  in  Panels  (a)  and  (b)  of  Figure  11.  As  can  be  seen  the  initial  Maxwellian  dis¬ 
tribution  function  of  the  electrons  remains  unchanged,  although  the  densities  (counts) 
do  change  with  time.  Panel  (c)  of  Figure  1 1  shows  the  density  profile  of  the  elec¬ 
trons  at  t  =  600o)pe-1 .  Note  that  the  fluctuations  are  as  much  as  50X.  The  fact 
that  the  temperature  and  the  distribution  function  of  the  electrons  remain  unchanged 
throughout  the  run.  further  substantiates  the  suggestion  made  earlier  that  the  elec¬ 
trons  move  in  a  coherent  manner  in  the  potential  field  of  the  waves.  Therefore,  as 
the  waves  grow  the  energy  of  the  electrons  also  increases,  and  as  waves  get  damped 
their  energy  decreases.  This  behavior  of  the  electrons  suggests  that  they  may  be 
approximated  as  a  charge  neutralizing  fluid  when  considering  the  ion/ion  mode.  Such 
an  approximation  is  desirable,  because  in  the  simulations  one  would  not  have  to  be 
concerned  with  the  electron  dynamics,  and  therefore  realistic  mass  ratios  between 
the  electrons  and  the  ions  can  be  used. 

3.3  Simulations  with  Boltzman  electrons 

To  approximate  the  electrons  as  a  charge  neutralizing  fluid,  use  can  be  made  of 
the  Boltzman's  relation 

ne(x)  =  n0  exp(e$(x)/Te)  (1) 


which  gives  the  electron  density  at  any  position  ne(x)  as  a  function  of  the  electro¬ 
static  potential  <t>(x)  at  that  position,  the  background  (or  unperturbed)  density  n0.  and 
the  electron  temperature  Te(where  the  Boltzman  constant  k  is  absorbed  in  Te).  The 
use  of  equation  (1)  is  expected  to  be  quite  satisfactory  when  the  density  fluctuations 
are  small  relative  to  n0.  and  therefore  can  be  thought  of  as  perturbations.  On  the 
other  hand  once  the  wave  amplitudes  become  sufficiently  large,  such  that  the  density 
fluctuations  are  a  large  percentage  of  n0.  it  is  not  obvious  if  and  why  relation  (1) 
should  hold.  In  fact  if  the  ion  density  (np  fluctuations  become  large  one  might  think 
that  a  better  approximation  for  the  electrons  may  be  to  replace  n0  in  equation  (1) 
with  n^x).  To  see  which  one  of  these  methods  better  approximates  the  actual  den¬ 
sity  of  the  electrons,  the  expressions  n0  exp(e$/Te)  and  nj  exp(e<f>/Te)  have  been 
computed  at  t  =  SOOcUpg-1  and  shown  in  panels  (d)  and  (e)  of  Figure  1 1  respectively. 
Note  that  in  our  case,  n^Cx)  is  the  sum  of  the  proton  and  lithium  densities  at 
position  x.  A  comparison  of  panels  (d)  and  (e)  with  (c)  clearly  indicates  that  even  in 
the  presence  of  large  amplitude  waves,  the  use  of  n0  provides  a  better  approximation 
to  the  density  of  the  electrons.  Note,  however,  that  even  in  this  case  notable  differ¬ 
ences  between  Panels  (c)  and  (d)  exist,  suggesting  that  the  Boltzman  approximation  is 
not  a  perfect  one. 

In  order  to  evalutate  the  usefulness  of  the  Boltzman  approximation  in  studying 
the  nonlinear  evolution  of  the  ion/ion  mode,  a  test  simulation  run  was  performed. 

This  run  was  similar  in  all  aspects  (i.e.  plasma  parameters,  mass  ratios,  etc.)  to 
the  full  particle  run  except  for  the  fact  that  the  electrons  were  approximated  by 
equation  (1).  The  full  nonlinear  Poisson's  equation  was  solved  for  using  the  Newton 
itteration  method  [Hockney  and  Eastwood,  1981].  The  phase  space  density  plots  of 
the  protons  and  the  lithium  ions  at  t  =  SOOcppg"1  from  the  test  run  are  shown  in 
panels  (a)  and  (b)  of  Figure  12  respectively.  A  comparison  between  this  figure  and 
panels  (c)  in  Figures  9  and  10  shows  that  aside  from  minor  details,  the  results  are 
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identical.  This  clearly  demonstrates  the  validity  and  the  usefulness  of  the  Boltzman 
approximation  in  describing  the  nonlinear  physics  of  the  ion/ion  instability. 

As  mentioned  earlier,  the  fluid  description  of  the  electrons  allows  one  to  use 
plasma  parameters  that  could  correspond  to  the  observation.  In  the  remainder  of 
this  section  the  results  of  a  simulation  run  with  realistic  plasma  parameters  are 
presented.  The  length  of  the  simulation  box  (L)  is  64  proton  Debeye  lengths  (Xp)  and  i; 

consists  of  64  cells.  Each  of  the  ion  species  are  represented  by  16,000  particles  and 
the  time  step  At  =  0.2o)pj-1  where  <opj  is  the  proton  gyrofrequency.  Figure  13 

shows  the  growth  rates  (20  of  the  ion/ion  instability  as  a  function  of  the  mode 
LK 

number  (-^-)  for  the  parameters  shown  in  the  figure.  These  parameters  are  simi¬ 
lar  to  those  used  by  Gurnett  et  al.  (1986a,b)  and  are  used  in  the  simulation  run.  As 
can  be  seen,  the  growth  rates  peak  at  mode  number  4  and  thus  one  would  expect  this 

mode  to  be  dominant  in  the  simulation.  Note,  however,  that  the  growth  rates  of 
modes  3  and  5  are  also  appreciable  and  thus  these  modes  may  also  show  up  in  the 
simulations.  Panel  (a)  of  Figure  14  shows  the  electric  field  of  the  waves  as  a 
function  of  x  at  t  =  606)pi-1  during  the  simulation.  It  is  evident  from  this  panel  that 
the  wave  amplitude  varies  with  x  or  in  other  words  is  being  modulated.  To  deter¬ 
mine  which  modes  have  been  excited,  a  Fourier  analysis  of  the  electric  field  squared 
has  been  performed  and  is  shown  in  panel  (b)  of  Figure  14.  As  can  be  seen  the  peak  in 
the  spectrum  corresponds  to  mode  4  with  modes  3  and  5  having  the  second  and  the 
third  largest  power  repectively.  This  result  is  in  good  agreement  with  the  linear 
theory.  Note  that  although  the  growth  rates  of  modes  3  and  5  are  comparable,  one 
would  still  expect  mode  3  to  have  larger  power  due  to  the  fact  that  the  initial  noise 
level  in  the  system  drops  with  increasing  K. 

Figure  15  shows  the  phase  space  density  plot  of  the  protons  at  three  separate 
times.  It  is  evident  in  panel  (a)  of  Figure  15  that  initially  the  protons  have  a  drift 
speed  Vd  =  2.2vp  which  is  smaller  than  the  beam  speed  of  the  protons  in  the  full 
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particle  simulation.  The  phase  velocity  of  the  unstable  waves  =  0.46vp  which 
falls  within  the  distribution  function  of  the  protons.  Thus,  one  would  expect  the 
trapping  of  the  protons  to  take  place  on  a  faster  time  scale  as  compared  to  the  full 
particle  simulation.  It  can  be  seen  in  panel  (b)  of  Figure  15.  that  four  trapping 
vortices  have  been  formed  as  expected  from  the  Fourier  analysis  of  the  waves.  Note, 
however,  that  the  two  vortices  on  the  right  hand  side  of  the  box  are  more  pronounced 
due  to  the  fact  that  the  electric  field  of  the  wave  is  larger  on  that  side  of  the  box. 
Panel  (c)  of  Figure  15  shows  the  protons  at  later  time  when  the  vortices  have  begun 
to  collapse  leading  to  the  scattering  of  the  protons  in  phase  space.  Note  that  the 
vortex  on  the  right  hand  side  corresponds  to  the  vortex  between  32  <  x  <  48  Xp 
shown  in  panel  (b).  Similarly  the  remnants  of  the  vortex  between  48  <  x  <  64  Xp 
in  panel  (b)  can  be  seen  at  0  <  x  <  1 6  Xp  in  panel  (c).  The  observed  speed  of  the 
vortices  is  near  the  phase  velocity  of  the  waves  as  one  would  expect. 

The  phase  space  density  plots  of  the  lithium  ions  are  shown  in  Figure  16  for  the 
same  times  as  those  in  the  previous  figure.  Imci3lly  these  ions  form  a  cold  non- 
drifting  Maxwellian  (panel  a).  As  the  waves  grow,  the  lithium  ions  are  acclerated  by 
their  electric  field,  leading  to  the  observed  fluctuations  in  panels  (b.c).  Note 
however,  that  these  ions  have  not  become  trapped  in  the  potential  well  of  the  waves 
in  contrast  to  the  protons  (see  panel  (b)  of  Fig.  15).  As  mentioned  earlier  this 
behavior  is  expected  since  the  trapping  time  for  the  lithium  ions  is  longer  than  that 
of  the  protons.  The  fact  that  the  phase  velocity  of  the  growing  waves  falls  within 
the  distribution  function  of  the  protons  results  in  their  rap;d  trapping,  and  satura¬ 
tion  of  the  instability  before  the  lithium  ions  get  trapped.  Despite  this,  however  it 
is  evident  that  some  of  the  lithium  ions  are  accelerated  to  relatively  large  veloci¬ 
ties  by  the  electric  field  of  the  waves. 

These  results  clearly  indicate  that  the  saturation  mechanism  of  the  ion/ion  mode 
remains  to  be  ion  trapping  when  realistic  plasma  parameters  are  used.  Nonetheless. 


differences  between  the  two  cases  also  exist,  in  that  the  effect  of  the  waves  on  the 
heavier  ions  is  not  as  pronounced  when  the  observed  plasma  parameters  are  used. 
Simulations  with  barium  ions  have  shown  similar  results  in  that  the  instability 
saturates  via  proton  trapping  and  the  much  heavier  barium  ions  are  not  strongly 

scattered  in  phase  space.  A  point  that  should  also  be  noted  is  that  the  observed 
heating  and  acceleration  of  the  ions  is  along  the  K  vector  of  the  ion/ion  wave  which 
is  oblique  to  the  solar  wind  flow  velocity.  Implications  of  these  results  in  regards 
to  the  role  of  the  ion/ion  mode  in  providing  coupling  between  flowing  plasmas  and 
the  new  born  ions  is  discussed  in  the  next  section. 


4.  SUMMARY  AND  DISCUSSION 


4.1  Summaru 

In  this  paper  the  nature  and  the  nonlinear  evolution  of  the  electrostatic  waves 
observed  during  the  AMPTE  solar  wind  releases  was  investigated.  By  using  the  full 
electromagnetic  dispersion  relation  the  properties  of  the  etectron/ion  and  the  ion/ion 
cross-field  instabilities  were  studied  and  their  growth  rates  were  calculated  for  the 
plasma  parameters  representative  of  those  during  the  AMPTE  releases.  It  was  shown 
that  these  growth  rates  are  an  order  of  magnitude  or  more  smaller  than  those  of  the 
electron/ion-acoustic  and  the  ion/ ion -acoustic  instabilities.  The  linear  analysis  of 
the  cross-field  instabilities  also  show  that  the  wavelength  of  the  maximum  growing 
wave  is  about  100  km  which  is  comparable  to  the  size  of  the  cloud  and  therefore  not 
likely  to  occur.  Similarly  the  frequencies  of  the  cross-field  modes  were  shown  to 
be  much  smaller  than  the  observed  frequencies.  By  using  the  magnitude  and  the 
direction  of  the  ambient  magnetic  field  measured  during  the  September  1 1  lithium 
release,  and  also  the  wave  electric  field  spectrums,  it  was  shown  that  neither  the 
frequency  nor  the  intensity  of  the  observed  waves  correlate  in  any  clear  manner  with 
the  changes  in  the  magnetic  field.  Since  one  would  expect  such  a  correlation  to  exist 
had  the  waves  been  generated  by  the  cross-field  instabilities,  it  was  clear  that  the 
observations  do  not  support  the  suggestion  by  Papadopoulos  et  al.  (1987)  that  the 
electrostatic  waves  were  generated  by  the  cross-field  instabilities.  In  view  of  both 
theoretical  and  observational  evidence  it  was  concluded  that  the  electrostatic  waves 
associated  with  the  AMPTE  solar  wind  releases  were  generated  by  the  ion-acoustic 
type  instabilities. 

Because  among  the  two  ion/acoustic  type  instabilities  the  ion/ion  mode  can  lead 
to  heating  and  acceleration  of  the  ions,  the  nonlinear  evolution  of  this  instabilitiy 
was  investigated  via  numerical  simulations.  By  using  a  full  particle  electrostatic 


code  it  was  illustrated  that  the  saturation  mechanism  of  the  ion/ion-acoustic 


instability  is  ion  trapping,  while  the  electrons  behave  as  a  charge  neutralizing  fluid. 
Because  of  this  characteristic  of  the  ion/ion  instability,  it  was  possible  to  conduct  a 
simulation  study  in  which  the  electrons  were  treated  as  a  Boltzman  fluid,  thus 
allowing  for  realistic  plasma  parameters  to  be  used.  The  results  of  this  simulation 
study  confirmed  the  conclusions  reached  by  the  full  particle  code  that  the  saturation 
mechanism  of  the  ion/ion  mode  is  ion  trapping  and  that  this  process  could  take  place 
using  realistic  prarmeters.  In  spite  of  these  similarities,  however,  differences 
between  the  two  simulation  studies  also  exist.  A  major  difference  being  that  unlike 
the  full  particle  simulations  in  which  both  the  protons  and  the  lithium  ions  became 
trapped,  in  the  Boltzman  electron  simulations  only  the  protons  became  trapped.  This 
contrast  was  attributed  to  the  fact  that  in  the  former  case  the  proton  drift  speed 
was  much  larger  than  the  wave  phase  speed  indicating  that  the  protons  would  first 
have  to  be  decelerated,  and  then  trapped;  thereby  increasing  their  effective  trapping 
period.  This  increase  in  the  trapping  period  of  the  protons  then  made  it  possible  for 
both  of  the  ion  species  to  interact  strongly  with  the  waves.  On  the  other  hand, 
during  the  fluid  electron  simulations  the  phase  velocity  of  the  unstable  waves  fell 
within  the  proton  distribution  function,  resulting  in  their  more  rapid  trapping. 

4.2  Discussion 

The  results  presented  here  clearly  show  that  the  ion/ion-acoustic  instability 
can  to  some  degree  result  in  transfer  of  energy  and  momentum  from  ions  of  a  flow¬ 
ing  plasma  to  the  newly  created  ions.  A  question  that  arises  then  is  how  effective 
this  mechanism  can  be  as  compared  to  the  other  modes  that  may  also  be  generated  due 
to  the  collective  plasma  processes.  Obviously  .the  answer  may  to  some  extent  depend 

'  on  the  situation  under  investigation.  For  example  it  was  suggested  by  Winske  et  al. 

! 

|  (1984,  1985)  that  during  the  AMPTE  releases  in  the  solar  wind  low  frequency  elec- 
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tromagnetic  waves  could  be  generated  by  the  ion-beam  instabilities.  However, 
because  of  their  long  wave  length  relative  to  the  cloud  size  these  waves  could  not 
have  been  generated  in  the  initial  stages  of  the  cloud  expansion,  and  thus,  the  ion/ 
ion-acoustic  mode  would  play  a  dominant  role.  On  the  other  hand  at  times  or  places 
(e.g.  comets)  where  the  electromagnetic  waves  can  grow  to  large  amplitudes  they 
become  very  effective  in  transfering  energy  and  momentum  between  the  ion  species 
(see  e.g.  Winske  et  al.t  1984  and  Omidi  and  Winske.  1987).  Whether  in  such  situa¬ 
tions  the  ion/ion-acoustic  instability  can  play  a  dominating  role  on  ion  dynamics  is  a 
subject  of  further  investigation.  However,  the  results  presented  in  Figure  16  do  not 
suggest  so.  For  example,  using  the  plasma  parameters  of  these  simulations  it  can  be 
shown  that  the  lithium  thermal  speed  corresponds  to  vL  =  1.84  km/s.  From  panel  (c) 
of  Figure  16  it  is  evident  that  some  lithium  ions  have  been  accelerated  to  velocities 
as  high  as  40  km/s  while  the  majority  of  the  ions  have  velocities  around  and  below 
20  km/s.  These  speeds,  however,  are  less  than  10*  of  the  solar  wind  flow  velocity 
and  therefore  not  considerable  in  terms  of  what  is  required  to  assimilate  the  ions 
into  the  solar  wind.  This  is  in  contrast  to  the  pick-up  by  the  elctromagnetic  waves 
which  can  result  in  average  flow  velocities  for  picked-up  ions  that  are  50*  of  the 
solar  wind  speed  (e.g.  Omidi  and  Winske,  1987).  Nonetheless  the  ion/ion-acoustic 
waves  may  also  play  an  important  role  in  the  pick-up  process.  To  further  investigate 
this  process  it  is  best  to  use  simulation  models  with  non-periodic  boundary  condi¬ 
tions  so  that  a  constant  refurbishing  of  the  solar  wind  protons  can  take  place.  The 
advantage  of  such  a  model  which  is  more  realistic  is  that  the  instability  may  not 
become  saturated  as  early  as  it  does  in  a  periodic  system.  Observations  of  the 
electrostatic  waves  during  the  AMPTE  releases  in  the  solar  wind  as  well  as  at  comet 
Giacobini-Zinner  (Scarf  et  al.,  1986)  are  good  encouragements  for  conducting  such 


studies  in  the  future. 
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Figure  5 


One  minute  measurements  by  the  magnetometer  and  the  plasma  wave 
instrument  during  the  September  1 1  lithium  release  in  the  solar  wind 
are  shown.  The  top  three  panels  show  the  direction  (azimuthal  and  lati¬ 
tudinal  angles  in  GSE  coordinates)  and  the  magnitude  of  the  magnetic 
field  while  the  bottom  panel  shows  the  electric  field  spectogram.  The 
onset  of  i.he  intense  electrostatic  waves  is  at  about  07:25:30  UT. 

The  growth  rates  of  the  ion/ion-acoustic  and  the  electron/ion-acoustic 
instabilities  maximized  over  the  wave  number  are  shown  as  a  function  of 
the  angle  (e)  between  the  K  vector  and  the  solar  wind  velocity.  The  solid 
line  represents  the  growth  rates  in  the  absence  of  cold  electrons,  while 
the  broken  line  corresponds  to  the  case  where  photoelectrons  are  present. 
For  the  plasma  parameters  shown  the  ion/ion-acoustic  instability  occurs 
at  0  >  70*. 

The  relative  orientations  of  the  solar  wind  velocity  (Vsw).  the  ambient 
magnetic  field  (B0)  and  the  K  vector  of  the  waves  in  the  (p,T\,C)  coordi¬ 
nate  system  are  shown. 

The  growth  rates  of  the  ion/ion  cross-field  instability  (<t>  *  88*)  and  the 
electron/ion  cross-field  instability  maximized  over  <  are  shown  as  a 
function  of  $  at  four  values  of  electron  0e.  In  this  figure  e  =  0°  and 
thus  K  lies  in  the  (p.O  plane.  Note  that  as  0e  increases  the  growth 
rates  become  smaller  due  to  electromagnetic  effects.  is  the  proton 
gyrofrequency. 

The  growth  rates  of  the  ion/ion  cross- field  instability  (panel  a)  and  the 
electron/ion  cross-field  instability  (panel  b)  maximized  over  1C  are  plot¬ 
ted  as  a  function  of  0  for  a  number  of  electron  0e.  Not©  that  the  growth 
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rates  of  the  ion/ion  mode  peak  at  larger  values  of  e  as  0e  increases 
while  the  maximum  growth  of  the  electron/ion  mode  occurs  at  e  =  0*  for 
all  values  of  0e. 

Figure  6  The  real  (solid)  and  the  imaginary  (broken)  parts  of  the  frequency  corres¬ 
ponding  to  the  ion/ion  cross-field  (panel  a)  and  the  electron/ion  cross¬ 
field  (panel  b)  instabilities  are  plotted  as  a  function  of  K  for  0e  =  °-5 
and  values  of  9  and  ♦  at  which  maximum  growth  occurs,  pj  is  the  proton 
gyroradius. 

Figure  7  In  this  figure  eight  electric  field  power  spectrums  covering  the  period 
07:25:30.2  -  07:25:38  are  shown.  Note  that  despite  the  large  changes  in 
the  magnetic  field  strength  during  this  period  the  frequency  at  peak 

power  remains  relatively  fixed  (f  >  100  Hz),  and  at  values  much  larger 
than  those  predicted  by  the  cross-field  instabilities. 

Figure  8  The  energy  history  of  the  electrostatic  waves,  electrons,  protons,  and  the 
lithium  ions  are  shown  from  top  to  bottom  panels  respectively.  Note 
that  the  energy  history  of  the  electrons  is  similar  to  that  of  the  waves. 

Figure  9  The  phase  space  density  plots  of  the  protons  are  shown  at  three  different 
times  during  the  simulations.  Initially  (panel  a)  the  protons  form  a 
beam.  At  later  times  they  get  decelerated  (b)  and  trapped  (c)  in  the 
electric  potential  of  the  waves  which  eventually  leads  to  their 
scattering  in  phase  space. 

Figure  10  The  phase  space  density  plots  of  the  lithium  ions  are  shown  at  three 

different  times  during  the  simulations.  As  can  be  seen  the  lithium  ions 
also  interact  strongly  with  the  excited  waves  leading  to  their  acceler¬ 
ation  and  eventual  scattering  in  phase  space. 

Figure  1 1  The  velocity  distribution  function  of  the  electrons  at  t  =  0  and  600o)pe_1 
are  shown.  As  can  be  seen  the  electrons  do  not  experience  any  heating 
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Figure  12 


Figure  1 3 


Figure  1 4 


Figure  1  5 


Figure  16 


due  to  the  wave-particle  interactions.  Panels  (c).  (d).  and  (e)  show  the 
electron  density.  n0  exp(e$/Te),  and  n^  exp(e$/Te)  as  a  function  of  x  at  t 
=  600o)pe-1.  Note  that  the  electron  density  profile  is  best  approximated 
by  n0exp(e$/Te). 

The  phase  space  density  plots  of  the  protons  and  the  lithium  ions  obtained 
from  a  simulation  code  with  Boltzman  electrons  are  shown  in  panels  (a) 
and  (b)  respectively.  A  comparison  between  these  panels  and  panel  (c)  in 
Figures  9  and  10  clearly  demonstrates  the  validity  of  the  Boltzman 
approximation  for  the  electrons. 

The  growth  rates  75  of  the  ion/ion-acoustic  mode  for  the  plasma 

parameters  shown  in  the  figure  are  plotted  as  a  function  of  mode  number 
KL 

( — ).  Note  that  the  growth  rates  peak  at  mode  number  4. 

2tt 

The  electric  field  profile  and  its  power  spectrum  at  t  =  60copj-1  during 
the  simulation  are  shown  in  panels  (a)  and  (b)  respectively.  Note  that  in 
agreement  with  linear  theory  the  largest  power  occurs  at  mode  number  4. 
The  phase  space  density  plots  of  the  protons  are  shown  at  three  different 
times  during  the  simulations.  As  can  be  seen,  similar  to  the  full  parti¬ 
cle  case  the  protons  become  trapped  and  later  scatter  in  phase  space  due 
to  interaction  with  the  excited  waves. 

The  phase  space  density  plots  of  the  lithium  ions  are  shown  at  three 
different  times  during  the  simulations.  Note  that  while  these  ions 
undergo  some  interaction  with  the  waves  they  do  not  become  trapped  in 
the  potential  field  of  the  waves. 
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